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Hermitian, symmetric and symplectic 

random ensembles: PDEs for the 

distribution of the spectrum 

By M. Abler and P. van Moerbeke* 

o 

psj ■ Abstract 

^ ■ Given the Hermitian, symmetric and symplectic ensembles, it is shown 

•^ . that the probability that the spectrum belongs to one or several intervals sat- 

isfies a nonlinear PDE. This is done for the three classical ensembles: Gaussian, 
Laguerre and Jacobi. For the Hermitian ensemble, the PDE (in the boundary 
points of the intervals) is related to the Toda lattice and the KP equation, 
^ . whereas for the symmetric and symplectic ensembles the PDE is an inductive 

equation, related to the so-called Pfaff-KP equation and the Pfaff lattice. The 
method consists of inserting time-variables in the integral and showing that 
Q^ ■ this integral satisfies integrable lattice equations and Virasoro constraints. 
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4. Proof of Theorems 0.1, 0.2, 0.3 

4.1. /9 = 2,1 

4.2. /3 = 4, using duality 

4.3. Reduction to Chazy and Painleve equations (/? = 2) 

5. Appendix. Self-similarity proof of the Virasoro constraints (Theorem 1.1) 

0. Introduction 

Consider weights of the form p{z)dz := e~'^^^'dz on an interval F = 
[A, B] C M, with rational logarithmic derivative and subjected to the following 
boundary conditions: 

(0.0.1) -^ = V'=^- = ^^^, \\^J{z)p{z)z^ = for all k > 0, 

together with a disjoint union of intervals, 

r 

(0.0.2) E = \J[c2^-l,C2i]CFCR. 

1 

The data (0.0.1) and (0.0.2) define an algebra of differential operators 

(0.0.3) ^, = ^cf+V(Q)^. 

Let Tin, Sn and 7^ denote the Hermitian (M = M ), symmetric (M = M ) 
and "symplectic" ensembles (M = M^, M = JMJ^^), respectively. Tra- 
ditionally, the latter is called the "symplectic ensemble," although the ma- 
trices involved are not symplectic! These conditions guarantee the reality of 
the spectrum of M. Then, 7in{E), Sn{E) and Tn{E) denote the subsets of 
Tin, Sn and 7^ with spectrum in the subset E C F CM.. The aim of this paper 
is to find PDEs for the probabilities 

(0.0.4) 

Pn{E) : = Pni all spectral points of M e E) 

^ Ih„(e), s„{e) or %,{E) g'*"" dM 

In„{F), Sn(F) or r„(F) e"**" ^'■^^dM 

= /..|A„(.)|^nL.e-(-)d..' /5 = 2, M respectively, 

for the Gaussian, Laguerre and Jacobi weights. The probabilities involve pa- 
rameters f3,a,b (see (0.1.1), (0.2.1) and (0.3.2)) and 



K4 



2 



^■-((r-(rMw::;;;M. 
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The method used to obtain these PDEs involves inserting time-parameters 
into the integrals, appearing in (0.0.4) and to notice that the integrals obtained 
satisfy 

• Virasoro constraints: linear PDEs in t and the boundary points of E, and 

• integrable hierarchies: 



ensemble 



Hermitian 

symmetric 
symplectic 



/3 lattice 






Toda 
Pfaff 
Pfaff 



As a consequence of a duality (explained in Theorem 1.1) between /3- Virasoro 
generators under the map (3 i-^ 4//3, the PDEs obtained have a remarkable 
property: the coefficients Q and Qi in the PDEs are functions of the variables 
n, /?, o, 6, and have the invariance property under the map 

h 



n 



to be precise, 
(0.0.5) 



Q^{-2n,f3,■ 



-2n, a 



'2' 2' 



:, h 



2' 



13=1 



Qi{n,(3,a,b)\p__ 



■A ■ 



Important remark. For (3 = 2, the probabilities satisfy PDEs in the bound- 
ary points of E, whereas in the case (3 = 1,4, the equations are inductive. 
Namely, for /3 = 1 (resp. /3 = 4), the probabilities Pn+2 (resp. Pn+i) are given 
in terms of Pn~2 (resp. Pn-i) and a differential operator acting on i-*„. 

0.1. Hermitian, symmetric and symplectic Gaussian ensembles. Given the 



disjoint union E C M. and the weight e ^^ 
on the form 



B, 



2r 

E 



cf+i 



the differential operators Bk take 



d 



da 



Also, define the invariant polynomials (in the sense of (0.0.5)) 



Q 



I2h^n [n + l-- 



Q2=4(l + <5f,J6(^2n + 5f,4(l-|)) 



and 



Qi 



'lA 



/3' 



Theorem 0.1. The following probabilities for {(3 = 2, 1,4) 



(0.1.1) 



Pn{E) 
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satisfy the PDE's {F := F„ = logP^): 

(0.1.2) 

^0^1 n- ^ 'n+ ^ \ . , . , I 2 when n is even and = 1 

Ot aQ ^ — - — 1 with index < . , . , ., i n a 

V n / I when n is arbitrary and p = 4 

BU + iQ2 + GBliF)Bl^ + 4Qi(3^2 _ 4^_^g^ ^ g^^)^ p 

0.2. Hermitian, symmetric and symplectic Laguerre ensembles. Given the 
disjoint union E C M^ and the weight z"e~^^, the B^ take on the form 

Also define the polynomials, again respecting the duality (0.0.5), 

3 

-n(n- l)(n + 2a)(n + 2a+l), for /3 = 1 

3 

-n(2n + l)(2n + a)(2n + a- 1), for /? = 4 

Q2 = [3/?n2-^ + 6an + 4(l-|)a + 3Uf_4 + (l-a2)(l-5f4), 



Qi = (^/3n^ + 2an + (1 - |)aj , Qo = 6(2 - <5f,4)('^ + ^), 



Theorem 0.2. The following probabilities 

..on p.p^ V|A„(z)|^nLi4e-''^''rf^^ 

satisfy the PDeA: {F := Fn = log Pn) 
(0.2.2) 



p2 " / 

^ n / 

bu - 2(5f,4 + 1)^:^1 

+ (Q2 + 6^2 ^F - 4(5^4 + l)B^^F)B^_, - 36l^{Qi - B^iF)B^i 
+ Q-i{3Bi - ABiB^i - 2Bi) + Qo(2^o^_i - Bo)]f 



^with the same convention on the indices n it 2 and n it 1, as in (0.1.2) 
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0.3. Hermitian, symmetric and sym,plectic Jacobi ensem,hles. In terms of 
E C [—1, 1] and the Jacobi weight (1 — z)"(l + z)^ , the differential operators 
Bk take on the form 

Setting 6o = a — b, 6i = a+6, we introduce the new variables, which themselves 
have the invariance property (0.0.5): 

r = ^{bl + {br + 2-[5f) s = ^boib, + 2-/3) 

4 
qn = ^{pn + bi + 2- p){pn + 61), 

and the following polynomials in g = qn,r,s, thus invariant under the map 
(0.0.5): 

(0.3.1) Q = —((s'^-qr + q'^f-Al.rs^-Aqs'^-As'^ + q'^r 
16 V 

Q^ = 3s^ - 3qr - 6r + 2q^ + 23q + 24, 

Q2 = 3qs^ + 9s^ - 4q^ r + 2qr + 4q^ + lOq^ , 

Q^ = 3gs2_^6s2-3gV + g3 + 4g2, 

^4 = Qs^ -'iqr-Qr + q^ + 22q + 2A = Qi + {<6s^ -q^ -q). 
Theorem 0.3. The following probabilities 

/^„ iA„(z)i'3 nLi(i - zkTii + zkfdzk 



(0.3.2) Pn{E) 



/[_,_y„ |A„(z)|/3 nLl(l - ^fc)"(l + ^fc)'rf^fc 

safe/y the PDE {F = F^ = logP„): 

/or /3 = 2: 
(0.3.3) 

' 2Sli + (g - r + 4)^-1 - (4^-1 F - s)fi_i + 3(?fi^ - 2qBQ + 8^0^' 

-4(g - 1)616-1 + (46_iF - s)Bi + 2(46_iF - s)BoB-i + 2qB2 ]F 
+AB\F (2B0F + sb^^f) = 
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for p= 1,4: 



(0.3.4) 



'P ,2P 2 



= (q + 1) ( 4gZ3li + 12(4S-iF - s)S'li + 2 (g + 12) (4e_iF - s)SoS-i 
+3q^Bl - 4 (g - 4) gBiB-i + g(4B_iF - s)Bi + 20gBoS^i + 2g^S2 ) F 



+ ( Q2B\ - sQiB-i + QiBo \ F + 48(B-iF)* - 4&s{B-iFf + 2Qi{B-iFf 

+ l2q^{Bv,Ff + lGq{2q-l){B\F){Bv,F) + 2A{q-l)q{BUFf 
+24^ 2B_iF -s][(q + 2)BqF + (g + ?,)B\F ] B-iF. 



0.4. ODEs, when E has one boundary point. Assume the set E consists 
of one boundary point c = x, besides the boundary of the full range. In that 
case the PDEs in the previous section lead to ODEs in x: 

(1) Gaussian (n x n) matrix ensemble (for the function (3 = 2, 1,4): 

fn{x) = -— logP„(max \i < x) 
ax i 



satisfies 

P?, 



/P 2P 2 ^ 



= /r + 6/;? + (4^(<4 - 2) + Q2J fn - 4^(<4 - 2)/„. 

(2) Laguerre ensemble (for j3 = 2,1,4): all eigenvalues Aj satisfy Aj > 
and 

fn{x) = x--logP„(max Aj < x) 
ax i 

satisfies (with / := fn{x)) 
(0.4.2) 

<4Q ( ^"'^I"^' - 1) - Uy - ^(<4 - 2) - Q,x - ?,51,q\ f 



x^f"'-i25l,-l)x'f" + 6x^f'^ 
t'^x 

1 



X f 4(<5f 4 + 1)/ - ^(<5f 4 - 2) - 2Qox -Q2 + 25f 4 + 1 ) /'. 
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(3) Jacobi ensemble: all eigenvalues Aj satisfy — 1 < Aj < 1 and 

fn{x) = (1 -x^)-— logP„(maxAi < x) 
ax i 

satisfies (with / := fn{x))'- 

for f5 = 2: 
(0.4.3) 2{x^ - iff" + 4(x2 - 1) [xf - 3/2) 

+ (iQxf - q{x^ - 1) - 2sx -r) f - f{4f -qx-s)=0 



for/?= 1,4: 

/P,2P 2 

(0.4.4) Q "^U"'^ -l 



4(q + l)(x^ - ly I -q(x' - 1)/'" + (12/ - gx - 3s)/" + 6q(g - 1)/' 
- (a.2 - 1)/' (^24/(g + 3)(2/ - s) + 8fq{5q - l)x ~ q{q + l){qx^ + 2sx + 8) + Qa 
+ / ( 48/^ + 48 f {qx + 2x~s)+2f (8q^x^ + 2qx^ - 12qsx - 2Asx + Q4) 
— g(g + l)x(3qx + sx — 2qx — 3g) + QiX — Qis 



For P = 2, fn{x) satisfies a third-order equation (of the so-called Chazy- 
type) with quadratic nonlinearity in /. Then /„ also satisfies an equation, 
which is second-order in / and quadratic in /", which after some rescaling can 
be put in a canonical form. Namely, 

Gauss gn{z) = b^^/'^ fn{zb^^/'^) + |nz, 

Laguerre, g„(z) = /„(z) + \{2n + a)z + \, 

Jacobi gn{z) := -^fnix)\x=2z-i - |^ + ^ 

satisfies the respective canonical equations of Cosgrove [11] and Cosgrove- 
Scoufis [12], 

• g"^ = -Ag'^ + 4:{zg' - gf + Aig' + ^2, (Painleve IV) 

. {zg"f = {zg' - g) (-Ag'^ + A^{zg' - g) + A^ + A^g' + A4, (Painleve V) 
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. {z{z - l)g"f = {zg' - g) (^V^ - ^g' {zg' - g) + A^ 

+Aig''^ + A^g' + ^4, (Painleve VI) 

with coefficients wliich will be determined in Section 4.3. Each of these equa- 
tions can be transformed into the standard Painleve equations. 

For (3 = 1 and 4, the inductive partial differential equations (0.1.2), (0.2.2) 
and (0.3.4) are new. For (3 = 2 and for general -E, they were first computed 
by Adler-Shiota-van Moerbeke [7], using the method of the present paper. 
For (3 = 2 and for E having one boundary point, the equations obtained here 
coincide with the ones first obtained by Tracy- Widom in [20] , who saw them to 
be Painleve IV and V for the Gaussian and Laguerre distribution respectively. 
In his Louvain doctoral dissertation, J. P. Semengue, together with L. Haine 
[14], were led to Painleve VI for the Jacobi ensemble, for (3 = 2 and E having 
one boundary point, upon subtracting the Tracy- Widom differential equation 
([20]) from the ones computed with the Adler-Shiota-van Moerbeke method 
([7]). As we shall see, the classification of Cosgrove [11] and Cosgrove-Scoufis 
[12], (A. 3) leads directly to these results. 



1. Beta-integrals 

1.1. Virasoro constraints for (3-integrals. Consider the data from (0.0.1) 
to (0.0.3) and the t-deformations of the integrals (0.0.4), for general /3 > 0: 
{t := (ti,t2,---) and c= (ci, C2, ..., C2r)) 

„ n 

(1.1.1) In{t,c;(3) := / \An{z)\f J] {e^^ '^'lp(.Zk)dzk) for n > 0. 

•^^" k=i 

The main statement of this section is Theorem 1.1, whose proof will be outlined 
in the next subsection. In Section 5 (Appendix), we give a less conceptual 
proof, which is based on the invariance of the integral (1.1.2) below, under the 
transformation Zi f-^ Zi + ef{zi)z^^ of the integration variables. The central 
charge (1.1.6) has already appeared in the work of Awata et al. [10]. 

Theorem 1.1 (Adler-van Moerbeke [2]). The multiple integrals 

„ n 

(1.1.2) In{t,c;(3) := / \An{z)\'' J] (e^^ '''h{zk)dzk) for n > 



fc=i 



and 



At " 

(1.1.3) In{t,c- -) := / |A„(z)|4/^ n (e^^*''^/>(^fe)d^fe) , /or n > 0, 
' •'^" k=i 
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with Iq = 1, satisfy respectively the following Virasoro constraints^ for all 
k > -1: 



i>0 ^ 



2 ' /3 



2 



"'" " •'^fc+i+l,n\^" 



(3t 2n 



)n/„(t,c;|)=0, 



in terms of the coefficients Oj, 6j of the rational function (—log/))' and i/ze end 
points Ci of the subset E, as in (0.0.1) to (0.0.3). For all n £ 1^, the ^Jl)j^(t,n) 

and Skn{t,n) form a Virasoro and a Heisenberg algebra respectively, inter- 
acting as follows: 



(1.1.5) 



^(2) /3^(2) 



ik-i)¥^l„+c 



12 



Jk,-e 






f3M 



me,n + ^'Hk + l)Sk,-i 



/3^(i) i3Aiy 

with central charge 
(1.1.6) 



c=l-6| (-1 






^y/2_ /^v'/'' 



and c 



1 1 

^"2 



Remark 1. The J^„'s are defined as follows: 

i+j— A; 

Componentwise, we have 



n '^'"k,n 



(1) 



/3 7(1)^ ^7(0) ..^ /3l(0) -^7(0) 



^(t, n) = ^4^' + njr and ^^ ^ = nJ^ = n5, 



Ok 



and hence 



^S^kl(t,n) 



I3\ (3j{2) 



'jP + 



n/3 + (fc + l)(l-^)) /54^) 



+ n((n-l)^ + l]jf, 



^When E equals the whole range F, then the B^'s are absent in the formulae (1.1.4). 
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(1.1.8) 


"4" = 


- el 4'-^"-' 




^4^> = 


,p 52 2 



^v ^u M ^ dt- 3^ ^^ ^ 

We put n explicitly in ^JI^j^(t,n) to indicate that the n component contains 
n explicitly, besides t. 

Remark 2. The Heisenberg and Virasoro generators satisfy the following 
duality properties: 

(1.1.9) Sjg(t,„) = ^J,g(-^.-|),„<,Z 

In (1.1.9), ^J^^ (— /3t/2, — 2n//?) means that the variable n, which appears in 
the rfi^ component, gets replaced by —2n/ (3 and t by —[it/2. 

1.2. Proof: 13-integrals as fixed points of vertex operators. The most 
transparent way to prove Theorem 1.1 is via vector vertex operators, for 
which the /3-integrals are fixed points. This is a technique which has been 
used by us already in [1]. Indeed, define the (vector) vertex operator X, for 

t = {ti,t2,...) e C°°, uGC: 

(1.2.1) X^(t,n) = A-ieSr*«"'e"'^^^'^^x(hl^), 

where xi^) •= (1) -2, ^^, • • •)• The vertex operator acts on vectors f{t) = 
(/o(t), /i(t), ...) of functions, as followsH 

(X^(t,n)/(t))^ = ei:r*«"" (|n|^)""Vn-i(t - /^K^]). 

For the sake of convenience, in this section we introduce the following vector 
Virasoro generators: ^T^'{t):= ( ^JJfc„(t, ?^))nez- 

Proposition 1.2. The multiplication operator z and the differential 
operators -^z^^^ with z G C*, acting on the vertex operator X^(f,z), have re- 
alizations as commutators, in terms of the Heisenberg and Virasoro generators 



q 2 3 

■^For a S C, define [a] := {a, ^, ^, ...) g C°°. The operator A is the shift matrix, with zeroes 
everywhere, except for I's just above the diagonal, i.e., {Av)n = v„^i. 
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(1.2.2) z''Xf3{t,z 

d 
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'j}\t),Xp{t,z) 



-/+'X,{t,z) 



^/^\t),Xp{t,z) 



Corollary 1.3. Given a weight p{z)dz on M satisfying (0.0.1), we have 

:i.2.3) 

d_ 
dz 



z'^'fiz)Xpit,z)piz) 



E («^ ^4'i(i) - h ^jS.+i(t)) , X^(t, z)p{z) 



i>0 



Proof. Using (1.2.2) in the last line, compute 
(1.2.4) -^^z''+^f{z)Xpit,z)piz) 



(pM. 
Vp(^) 



f{z)^ z^+^Xp{t,z)p{z)+p{z)^^ (z^+^f{z)Xp{t,z)) 



d 



E hz'^+'+^Xpit, z) p{z) + p{z)— E a./+*+^X^(t, z) 





oo 



dz 



L 



Eft.''4+m.X;3(t,^)p(^) + E^^'^jfi'X^li'^M 



L 



establishing (1.2.3). 



D 



Given the weight pE{u)du = p{u)lEiu)du, with p and E as before, and 
with Ie the indicator function of E, define the integrated vector vertex operator 



(1.2.5) Yp{t,pE) 
and the vector operator 

(1.2.6) Vk := Bk - Vk 

2r 



dup{u)Xp{t, u), 



d 



1 '^^^ j>0 

consisting of a c-dependent boundary part B]^ and a (i, n)-dependent Virasoro 
part Vfc. 

Proposition 1.4. T/ie following commutation relation holds: 

(1.2.7) [Pfc,Y^(t,ps)]=0. 
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Proof. Integrating both sides of (1.2.3) over E, one computes: 



(1.2.^ 



2r 



^k+1 



^dz—[z'^'f{z)^p{t,z)p{z) 



1 

^ OCi Je 

[Bk,Y(,{t,pE)]; 



while on the other hand 



(1.2.9) / dz 

IE 



E («^ ^^kl - b^ ^S^kl+l) ,^p{t,z)p{z 



j>0 



h''s'il+,), f dzp{z)Xf,{t,z) 



E («^ '^'kl - 

= [Vk,Yp{t,pE)]. 
Subtracting both expressions (1.2.8) and (1.2.9) yields, using (1.2.3), 

0=[Bk- Vk,Yp{t,pE)] = [Vk,Ypit,pE)] , 
concluding the proof of Proposition 1.4. 

Proposition 1.5. The column vector, 

I{t) := U^ |A„(z)|'^ n eT.r^^^lp{zk)dzk 



n>0 



is a fixed point for the vertex operator Y^(t, pe)'- 
(1.2.10) (Y^(t,ps)/)„ = /„, n>l. 

Proof. We have 
(1.2.11) 

„ n 

In{t) = / \An{z)fY[(e^T^^^lpj,(zk)dZk 



/ dupEi 



u]e 

n—l 



k=l 

E;"t,n'| |/3(n-l) 



D 



„ n—l 13 n—l 

J^'-' k=l " k=l ^ 
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/dnp£(n)eSr*>«>|/5(«-i) 



Y^(t,Pi,)/(t) 

It suffices to do the above argument for all tj > 0, enabling one to replace 
„ 1 '^ by |e^i '^ |. Then one continues the result for all tj G C. D 

Proof of Theorem 1.1. From Proposition 1.4 it follows that for n > 1, 

(1.2.12) = [Vk,{Yp{t,pE)T]I 

= VkYp{t,pETl-Yp{t,pETVkI. 

Taking the n component for n > 1 and fc > — 1, setting 

X/3(t,u) = eS^'^'e"^^"^^, 
and using (1.2.10), we have 

= {VkI-Yp{t,pETVkI)^ 

= {Vkl)n - J dupE{u)Xp{t- U){\ufr-\..J dupE{u)Xp{t- u){VkI\ 
= {Vkl)n- 

Indeed {T>kI)o = for /c > — 1, since /q = 1 and V^ involves Bkj^Jj., ^ J^ 
and J^ for A; > —1: 

Bk and ^ J^. are pure differentiations for /c > — 1; 

^ J^. is pure differentiation, except for k = —1; 

^ J_l appears with coefficient n(5, which vanishes for n = 0; 

J^ appears with coefficient n{{n — 1)^ + 1), vanishing for n = 0. 

The proof of the 2 formula in (1.1.4) follows immediately from the du- 
ality (1.1.9). D 

1.3. Examples. Example 1 {Gaussian [3- integrals) . The weight and the a^ 
and 6j, as in (0.0.1), are given by (setting 6 = 1 in (0.1.1)) 

p(^) = e-^(^) = e-^', V' = g/f = 2z, 

oo = 1, 6o = 0, 6i = 2, and all other Oj, hi = 0. 
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From Theorem 1.1, the integrals 

„ n 

(1.3.1) /„=/ \^n{z)f\{e-''"^^Z^'^''^dzk 

satisfy the Virasoro constraints 
(1.3.2) 

-B.In = - E cr^|-/„ = (- '^jgi + 2 '^JK,,.) In, fc = -1, 0, 1, . . . . 

Introducing the fohowing notation 

ai = {n- ^^)/3 + i + 1 - 6o = (n - ^^)/3 + i + 1, 

and upon setting F = login we find that the first three constraints have the 
following form: 

-so^ = i2|-|:«.|-)^-^. 

For later use, take linear combinations such that each expression contains 
the pure differentiation term dF/dti: 

(1.3.3) Pi = -l^_i, V2 = -^Bo, P3 = -^(6i + yfi-i 

which yields 



V2F = \- yiti^—\F cTi, 



dt2 2^^ 'du j 4 



V3F = \- > itjT; -7(^1} iti^ \ F - -aih. 
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Example 2 (Laguerre (3 -integrals). Here, the weight and the Cj and 6j, as 
in (0.0.1), are given by (again setting 6 = 1 in (0.2.1)) 

e-V' = ^«e-^ F' = 4 = i^^, 
/ z 

oo = 0, oi = 1, 6o = — O) bi = 1, and ah other Oj, hi = 0. 
Thus from (1.1.4), the integrals 

„ n 

(1.3.5) /„=/ |A„(z)|'3[]z^-^'=+E.=i*«4d^, 
satisfy the Virasoro constraints, for /c > — 1, 

(1.3.6) -B,I^ = - E C^ |-/n = (- ^Jgl,n - « ^jSl,n + ^tfi.n) ^- 

Introducing the following notation, as before, 

ai = {n- ^^)/3 + i + l-hQ = {n- ^-^)P + i + l + a, 
and upon setting F = Fn = log /„ , we see that the first three have the form: 

-fi-.^ = (j|-|:".|)f-f(-.-). 

Replacing the operators Bi by linear combinations Pj, we see that 

(1.3.7) Vi = -B.i 

V2 = -Bo-aiB-i 

2?3 = -Bi - (T2-B0 - aia2B^i 

yields expressions, each containing a pure derivative dF/dti 

(1.3.8) 

^ ^ dF ^^ dF n, 

^^^ = atr-g^*^a^-2(^^ + ")' 
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^ r. dF ^^ . ( d d d \^ n, 



^*3 Ki V dti dti+i dti+2 



2 [dtj +UtJ )' 

Example 3 (Jacobi (3-integral). The weight and the Cj and 6i, as in (0.0.1), 
are given by 

(\ -V (-, ^an , \b ^ri 9 a - b + {a + b)z 

Pab{z):=e ={l- z) [l + z) ,V =- = ^ > 

/ I- z^ 

ao = 1, ai = 0, 02 = —1, 6o = a — 6, 6i = a + 6, and all other a^, 6j = 0. 
The integrals 

„ n 

(1.3.9) / \^n{z)\^\{{l-ZkT{l + Zkfe^t.^^''^dzu 
■^E" k=i 

satisfy the Virasoro constraints (k > —1): 

(1.3.10) -Bkin = -flcl+H^-ch-^^In 

Introducing the following notation, 

i + 1 
ai = (n —)(3 + i + l + bi, 

and upon setting F = Fn = log/„, we see that the first four have the following 
form: 

(1.3.11) 



-^oi^ = C72— + 6o7^ + E^** 7^7 7^ +^7^^^ 

I 3t2 OTl ^ 3tj+2 OTj 2 9tf ' 



+^(^1 -^i^i-bi), 



" 2 V^tT/ "2' 



92 \ r,d£dF_ 

^^dhdr2dt'idh' 
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d , d , , ^ d ^ , d d 

i>l 



-"^^ = l'««i + ''»at^-<"^-'')a5 + g*W-8^' 






2. Matrix integrals and associated integrable systems 

2.1. Hermitian matrix integrals and the Toda lattice. Given a weight 
p{z) = e~ '^-^ defined as in (0.0.1), the inner-product 

(2.1.1) (/, g)t = I f{z)g{z)pt{z)dz, with pt := e^T ^^^'p(z), 

JE 

leads to a moment matrix 

(2.1.2) mn{t) = {Pijit))o<i,j<n = i{z\z^)t)o<i,j<n, 

which is a Hdnkel matrioO, thus symmetric. Hankel is tantamount to Arrioo = 
m-ooA . The semi-infinite moment matrix m.oo evolves in t according to the 
equations 



^^'^ -//■,. and thus ^"^^ - K^m I commuting 

dti. oth \ vector helds 



Another important ingredient is the factorization of ttt-oo into a lower- times 
an upper-triangular matrixB 

moo{t) = S{t)-^S{t)^-\ 

where S{t) is lower-triangular with nonzero diagonal elements. 

Theorem 2.1. The vector T{t) = {Tn{t))n>o, with 

1 /■ " 

(2.1.4) Tn{t) := detm„(t) = - / Al{z) IT pt{zk)dzk 

satisfies: 

(i) Virasoro constraints (1.1.4) for /? = 2, 

(2.1.5) [-E-rV(c.)|- +E (a, J^ - h J« +0 j r = 



"^Hankel means fiij depends on i + j only. 

^This factorization is possible for those t's for which Tri{t) := detmn(t) ^ for all n > 0. 
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(ii) the KP-hierarch 

of which the first equation reads: 



Tn°Tn = 0, 



dtj ^ [dt2) dtidts 



dt'i 



log T„ + 6 I ^g^ log Tn\ =0, 

A: = 0,l,2,... 



(iii) The standard Toda lattice; i.e., the tridiagonal matrix 

1/2 „ / xl/2 



(2.1.6) L{t) := S{t)KS{t)-^ 



(roT^\ _9_lng.Il fnsV 







(^) 



1/2 



7|-log^ 
ail ^ T2 



satisfies the commuting equations^ 



(2.1.7) 






:(^')s,^ 



(iv) Orthogonal polynomials: The n degree polynomials Pn{t; z) in z, 
depending on t £ C°°, orthonormal with respect to the t- dependent 
inner product (2.1.1) 

{pk{t;z),pi{t;z)) = 5ki 

are eigenvectors of L, i.e., {L{t)p{t; z))n = zpn{t; z), n > 0, and enjoy 
the following representations 



( 



Pn(t-z):={S{t)x{z)\ 



\J rn{t)Tn+\{t) 



:det 



mr. 



V Ain,0 



Tn+l{t) 



1 \ 

Z 



z" / 



Tn{t) Tn{t) 



"for the customary Hirota symbol p(dt)f o g := p{-^)f(t + y)g{t — y) 



The p/'s are the 



y=0 



elementary Schur polynomials e^i ' := y^.^„ Pi(ti, t2, . . .)z' and pf (9) := p^(^-, ^ ^-, • ■ ■)■ 
''Os means: take the skew-symmetric part of () in the decomposition "skew-symmetric" 
"lower- triangular . " 
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The functions qn{t; z) := z J^„ ^"J-J^' pt{u)du are "dual eigenvectors^^ 
of L, i.e., {L{t)q{t; z))n = zqn{t;z), n > 1, and have the following 
T-function representation: (see the remark at the end of this section) 



(2.1.8) qnit;z) := z f P^l^^p,{u)du = (s'^-\t)x{z-')) 

JM.n Z — U V / n 

= (s{t)mUt)x(z~^: 

^ ""71 /,\ 

rn{t) 

(v) Bilinear relations: for all n,m >0, and a,b G C°^, such that a — b 
t-t\ 



2-Ki 

dz 



(2.1.9) / Tn{t-[z'^])Trn+l{t' + [z-^])e^^''^~ Z 
J z=oo 

Jz=o ^iri 

In the case /? = 2, the Virasoro expressions take on a particularly elegant 
form, namely for n > 0, 

Jgi(*) = E ■ 4]lit) jS(i) : = 4'\t) + 2njf )(t) + n%, 

i+j=k 

4!iW = ji'\t) + n6ok, 
with! 

O 1 

(2.1.10) 4') = — + -(-k)t^,, 



. -f^ , dtidtj ^ , dtj A f^ , 

i+j=k -^ —i+]=k •> —i—j=k 



Statement (i) is already contained in Theorem 1.1, whereas the other state- 
ments can be found in [1], [2], and [5]. Notice that the standard Toda lattice 
is a reduction of the semi-infinite 2- Toda lattice, where Tn{t,s) = Tn{t — s). 
The 2- Toda lattice arises in the context of a factorization of a generic semi- 



infinite matrix m^{t, s), satisfying the simple equations ^°° = A'^ 



dmao A fc„on 9"i 



mr. 



—rrioo^^, whereas the standard Toda lattice is related to the same factoriza- 
tion of moo{t,s), but where moo{t,s) is Hankel (i.e., AttIoo = mooA.'^)- 



'The expression j). = for A; = 0. 
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Remark. The vectors p and q are eigenvectors of L. Indeed, remembering 

x{z) = (1, z, z'^, ...)"'', we have 

^X{z) = zx{z) and A^x{z'^) = zx{z~^) - zei, with ei = (1, 0, 0, ...)^. 
Therefore, p{z) = Sx{z) and q{z) = S^^^xiz^"^) are eigenvectors, in the sense 
Lp = SA.S^^Sx{z) = zSx{z) = zp, 
L^q = S'^-^K^S'^S'^-^xiz'^) 

= zS'^-^x{z'^) - zS'^-^ei = zq- zS"^~^ei. 
Then, using L = L^, one is lead to 

((L — zl)p)n = 0, for n > and ((L — zl)q)n = 0, for n > 1. 

2.2. Symmetric / symplectic matrix integrals and the P faff lattice. Consider 
an inner-product, with a skew-symmetric weight p(y,z), 
(2.2.1) 

{f,9)t = f{y)g{z)e>^T t^(y'+'') p{y, z)dy dz, with p{z,y) = -p{y,z). 



5 . . .y ; 



Then, since 

{f,9)t = -{g,f)t 

the (semi-infinite) moment matrix, depending on t = (ti,t2 

mnit) = {pij{t))o<i,j<n~l = {{y\z^)t)o<i,j<n~l 

is skew-symmetric and the semi-infinite matrix nioo evolves in t according to 
the commuting vector fields 

(2.2.2) -^ = la+kj + f^i,j+k, i.e., -^ = A^moo + rriooA^^ 

It is well known that the determinant of an odd skew-symmetric matrix equals 
0, whereas the determinant of an even skew-symmetric matrix is the square of 
a polynomial in the entries, the Pfaffian, with a sign specified below. So 

det(m2„_i(t)) = 

(detm2„(t))^/^ = pf{m2n{t)) = —{dxoAdxiA... Adx2n-iy^ 



n\ 



y^ jiij{t)dxi A dxj 



Define now the Pfaffian r-functions: 

(2.2.3) T2n{t):=pfm2n{t), 
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and the semi-infinite skew-symmetric matrix, everywhere, except for the 2x2 
blocks, along the diagonal: 



/ 



(2.2.4) 



J :-- 



1 






-1 








1 






-1 








1 






-1 



\ 



V 



with J2 = -/. 



/ 



Since rrioo is skew-symmetric, rrioo does not admit a Borel factorization in 
the standard sense, but moo admits a unique factorization, with the matrix J 
inserted (see [6]): 

m^{t) = Q-\t)JQ^-\t), 



(2.2.5) 


/ . 







Q{t) = 




Q2n,2n 
Q2n,2n 








* 


Q2n+2,2n+2 

Q2n+2,2n+2 



V 



G K. 



■■ J 



K is the group of lower-triangular invertible matrices of the form above, with 
Lie algebra 6 of matrices of precisely the same form. In this problem, the Lie 
algebra splitting of semi-infinite matrices is given by 



(2.2.6) gl{oo)=t®n' 



t = {lower-triangular matrices of the form (2.2.5)} 
n = sp(oo) = {a such that Jo' J = a}. 



with unique decomposition {a± refers to projection onto strictly upper- (strictly 
lower) triangular matrices, with all 2 x 2 diagonal blocks equal to zero) 



(2.2.7) 



(a)j + (a)n 



[{a. - J{a+)^J) + ^(ao - J(ao)^ J)) 



+ ( (a+ + J{a+) J) + -(ao + J(ao) ' J) ) . 
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Considering as a special skew-symmetric weight (2.2.1), 

(2.2.8) p{y,z) := 2D''6{y - z)p{y)p{z) ,with a = ^1, p{y) = e'^^^), 

the inner-product (2.2.1) becomeso (see [8]) 

{f,9)t = [ [ f{y)g{z)eT.^^^y'+-^'hD^5{y - z)p{y)p{z)dydz 



J{y)g{z)e^T *'(s''+^')e(y - z)p{y)p{z)dy dz, for a = -1 

{f.9]{y)e^^^''''^P{yfdy, for a = +1, 

and (see [16], [4]) 
(2.2.9) 

TTTTT / \^2n{z)\ n e'>^T''-lp{zk)dzk 
(2n)! jR2n ^J-_^ 

= r^ f eT'-(-^W+E*-^OdX, for a = -1, 
(2n)! Js2n 

-[ \/\n{z)tf{eT.r^^^^lp\z,)dzk 
n! Jr" ^J-^ 

1 /■ ^Tr(-mx)+Y.^UX^)dX, fora = +l. 



Setting 



n\ JT2n 

p{z) = p{z)Ie{z) for a = — 1 

p(2) = pi/2(^)j^(^)^ t ^ i/2 for a = -hi 



in the identities (2.2.9), we are led to the identities between integrals and 
Pfaffians, which are spelled out in Theorem 2.2: 



Theorem 2.2. The integrals In{t,c), 



n 

In = I \Aniz)fl[[e^T^^^lp(zk)d. 
^" k=i 



Zk 



^e{y) = sign(y), and {f,g} := f'g ~ fg'. Also notice that e' = 25(x). 
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nlpf ( / / y'z^iy - z)eSr ^'=(y''+^'^ p{y)p{z)dydz 

lE^ / 0<jj<n-l 

n even, for (3=1 
n!r2„(t/2,c), 



= n\Tn{t,c), 

W (^Jjy\y'}e^^''''^''p{y)dy 



0<i,j<2n-l 



n arbitrary, for (3 = A 



and the Tn{t,c)^s above satisfy the following equations: 
(i) The Virasoro constraintstj (1.1.4) for /? = 1,4, 



2r 



(2.2.10) I - E crV(Q)^ + E («•' "J^.n - b. ^J^^i J I In = 

1 '^'^« i>0 

(ii) The Pfaff-KP hierarchy: (see footnote 6) 

/ 1 a2 ^ 



(2.2.11) 



Pk+4{d) 



2dtidtk+zj 
of which the first equation reads 



Tn°Tn= Pk{d) r„+2 ° T„_2 

n even, fc = 0, 1, 2, ... 






(iii) The Pfaff lattice: T/ie time- dependent matrix 



12 



T"n-2Tn+2 



n even. 



(2.2.12) L(t) = Q{t)KQ{t)-^ 

satisfies the Hamiltonian commuting equations, given by the Adler- 
Kostant-Symes splitting theorem, applied to the splitting gl [oo) = t®n, 
as in (2.2.6) and (2.2.7), 

BT 

— = [-{L'% L], (Pfaff lattice) 

Oti 
(iv) Skew-orthogonal polynomials: The vector of time- dependent polyno- 
mials q{t; z) := {qn{t', z))n>o = Q{t)x{z) in z satisfy the eigenvalue 
problem 



(2.2.13) 



L{t)q{t,z) = zq{t,z) 



'^here the a^'s and bi's are defined in the usual way, in terms of p{z); namely, — — — -^^ 



E". 
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and enjoy the following representations: 

„ (f.^^ _ ^2nr-l/2 T2n{t-[z-^]) _ T2n+2{t) 

Q2n{t,Z) - Z ft,2n ^ /.x ) "-2n - ^ /,s 

T/iey are skew- orthogonal polynomials in z; i.e., 

{qi{t;z),qj{t;z))t = Jij. 

(v) The bilinear identities: For all n,m>0, the T2nS satisfy the following 
bilinear identity 



(2.2.14) / r2„(t - [.-^])r2„+2(t' + [z-i])eEr(*-*:y ,2n-2„.-2^ 
Jz=oo 2Tri 

+ i T2n+2{t + N)T2™(t' - [z])eT.7 i^'.-^^)^-\^^-^-^ ^, = Q. 

Jz=Q 2ttz 



Note that (2.2.10) is a consequence of Theorem 1.1, while items (ii) to 
(v) are shown in [4], [6]. (See [8] for the Pfaff lattice, viewed as a reduction 
of the 2-Toda lattice.) A semi-infinite matrix moo{t,s), satisfying !^°° = 

A^moo, ^§^ = —niooh^^, leads to the semi-infinite 2-Toda lattice. When the 
initial condition ?7ioo(0,0) is skew-symmetric, then moo{t,—t) remains skew- 
symmetric in time and r„(t) = {Tn{t, — t))^'^ = pfmn{t, —t) is a Pfaff lattice 
r-function. 



3. Expressing t-partials in terms of boundary-partials 

3.1. Gaussian and Laguerre ensembles. Given first-order linear operators 
Pi,2?2,^3 in c = (ci,...,C2r-) G M^'' and a function F(t,c), with t G C°^, 
satisfying the following partial differential equations in t and c: 

dF 
(3.1.1) VkF=--+ J2 7kjVj{F)+-fk + 6kti, fc= 1,2,3,..., 

'' -i<j<k 

with Vj (F) nonlinear differential operators in ti of which the first few are given 
here: 

(-. .,.,.^E;.|^.f..(f .(f)> -...^. 

In (3.1.1) and (3.1.2), /? > 0,7^^,7^,5^ are arbitrary parameters; also 62J = 
for J 7^ 2 and 52j = 1 for j = 2. The claim is that the equations (3.1.1) enable 
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one to express all partial derivatives, 



(3.1.3) 



dt\K..dt 



«fc 



, along C := {all tj = 0, c = (ci, ..., C2r) arbitrary}, 



uniquely in terms of polynomials in Dj-^...Dj^F{0, c). Indeed, the method con- 



sists oi expressmg — — 



in terms of T>kf 



t=o 



using (3.1.1). Second deriva- 



t=o 



tives are obtained by acting on V^F with 2?^, by noting that 2?£ commutes 
with all t-derivatives, by using the equation for ViF, and by setting in the end 
t = 0: 

V.VkF = Ve—+ Y. lkP,{V,{F)) 

^ -l<j<k 
d \ 

— h ^ IkiVj T^e{F), provided Vj{F) does not 



d_ 
dtk 

d^F 



'l<j<k 



-l<j<k 



dF 



contain nonlinear terms 

+ E ^ejVjiF) + 5eti 
-i<j<e 



+ lower- weight terms. 



dtkdti 
When the nonlinear term is present, it is taken care of as follows: 

_ fdF\^ dF^ dF 

dtidti 

f)F r) I clF 1 

higher derivatives are obtained in the same way. Explicit expressions for only 
a few partials, useful in the next subsection, will be given here: 



(3.1.4) 



dF 

d^F 



dtj 



ViF--fu 

(t>1 - 710^1 j F + 71071 - 6i, 



174 



M. ADLER AND P. VAN MOERBEKE 



d^F 



dtf 



dtf 
dF 

d^F 



Vf - 37io^f + 2710^1 ^ + 27io(<Ji - 7i7io 



Vt - 67ioPJ + ll7io^f - 6710^1 F - 6710(^1 - 71710) 



^2-^-72, 



9^2 



dF 



d^F 



Vl - 272o2?2 + /3721732I?? 

- ((271 + 710)721732/3 + 272 _i)Pi - 2721P3 ]F 

+ /3721 732(^1^)^ + /372i732(7i + 7io7i - h) 
+ 2(72173 + 72072 + 7172,-1), 



^3 - ^7321?? + f 732(271 + 710)^1) F - T:732(PiFf 



dhdts 



+ i7732(<Ji - 71710 - 7i) - 73, 



^1^3 - ^732^?? + /3732(7i + 2710)^? 



3/3 \ 

- -^710732(271 +710)^1 - 37i,_-iP2 - 3710^^3 j-P 

+ Y7io732(:Pii^)' -/3732(:Pii^)(I??i^) 

3 

+ 2 (^71073 + /37327io(7i + 71071 - ^1) + 271,-^172)- 



3.2. Jacobi ensemble. 

1. Prom the expressions (1.3.11), upon evaluating 13^iF\^^q , B^iF\ „ , 
BqFI^^q , one finds the following equations, both sides of which are evaluated 
at t = 0, 



1 
0-1 



-B^iF = ai--F + bon, 
oti 
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-BnF 



, d d 

Oti Ot2 



F + 



/3 



_d_ 
dti 



F + 



dF\' 



n 



((Ti - bi] 



From these expressions, one extracts 



dF_ 

dti 



t=o 






t=0 



dF 



t=o 



in terms of B^^F. 



2. From the expressions for B'^iF\^^q , BQB-iF\f.^Q , BiF\f.^Q, namely 



BiF 






(y\ 



BoB-iF 



1 



Bl^F 



0-1 
one extracts 






■dtidt2 dts dti dt{ 



2 \ dti dti dti 



2 ' 



^3 



a2 



^^^a^ + ^^^9tT9^-^Jr + ^J^]^' 



OF 

dt^o 



d^F 

t=o' dti 



d^F 



t=o 



' dtidt2 



t=o 



in terms of B^jF, using the previous extractions. 



3. From the expressions for B2F\^^q , BiB-iF\^^q , BqF\^^q , BoB'^iF\^^q, 
fiix^li=0' i^S'^sly, (where both sides are evaluated at t = 0) 

(3.2.1 

dt4 dti 



( d d d / d-' d'^ 9^ 



dt2 2\dtl dti 'dtidtz 



F 



+ P 



dF 



dF 
df2 



2 



dF dF \ 
dtlduj ' 



—BiB-iF 



B^oF 



_d d_ 9' ._£!_ 

dt4 dt2 dtidt2 ' dtidti 



d^ d^ 



+ /3 



d^FdF dF d^F 



dti dt2 dti dtidt2 



, d d 13 d'^ n d r. 9 

dti dt2 2 dtf dti dti 

2 



, dF dF 

'°dri+'''dr2 + 



\^ g-F dF. P d'F (dF\ 

2^'^'^dt,+2 dt?^ 2\dtl^\dti) 
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CTl oil V OTl OT2 



I ^ dF dF yr^ . , dF dF . (3 d'F (dF\ 

„4^ d f d d \ f d d ^ , d d , 



OF dF 



Of Ot v-^ . C'i' 



9F 9F 



at. 



) + ban — nti 



one extracts 
(3.2.2) 



d^F 



dtf 



t=o 



dF 



d^F 



' f)+2i 



d^F 



t=o 



' dhdt3 



d^F 






t=0 



t=o dt{dt2 
again in terms of S!i^, using all the previous extractions. 

3.3. Evaluating the matrix integrals on the full range. The denominators 
of the probabilities (0.0.4), for (3 = 1,4; namely: 



( t " 



j(/3) ., 



n 

^n{z)\^Xlzle-''^''dzk 



k=l 



f \Aniz)f flil - ZkTil + ZkfdZk 



can be evaluated, using Selberg's integral (see Mehta [16, p. 340]): 



/(/3) 



n— 1 



"• ^ "■ ' /J, r(/3/2 + l) 



n-l 



^-n(/3(n-l)+2a+2)/2 TT 

i=0 

2n{2a+2b+/3(n-l)+2)/2 



r(a + l+i/3/2)r((i + l)/3/2 + l) 
TiP/2 + 1) 



rt— 1 

n 



r(a + J/3/2 + l)r(b + J/3/2 + l)r((j + 1)13/2 + 1) 
r(/3/2 + l)r(a + 6 + (n + j - l)/3/2 + 2) 



Lemma 3.1. For future use, the following expressions 



' n(n— 1) 



6l^=^) := 



(n!) 



n2 



r(l) r(l) 
^n-2-'n+2 



(Gauss) 



(n-2)!(n + 2)! (jW) 



n(n- 1) (ra+2a) (n+2a+l) 

(Laguerre) 



L Qi 



(Jacobi) 
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6l^=^) 



(n!) 



1^2 



r{4) .(4) 
-'n-l-'n+l 



2n(2n+l) 
4P 



(Gauss) 



(n-l)!(n + l)! (jW)2 



2n(2n+l)(2n+a)(2n+a-l) 
P 

(Laguerre) 
-^ (Jacobi) 



satisfy the following functional dependence: 

In the expressions above, Q {already appearing in (0.3.1)), and a new expression 
Qq are expressible in terms of the variables q,r,s introduced in (0.3.1): 

48(n - l)n(2o + n)(2a + n + 1)(26 + n) (26 + n + 1) 

(26 + 2a + n + 1) (26 + 2a + n + 2) , for (/? = 1) 

96n (2n + 1) (a + 2 n - 1) (a + 2n) (6 + 2n - 1) 

(6 + 2 n) (6 + a + 2 n - 2) (6 + a + 2 n - 1) , for {P = 4) 



Q 



3 

16 



TE {{s^ - qr + q^f - 4(rs^ - Aqs^ - As^ 



\) 



an 



Qt 



48 (6 + a + n) (6 + a + n + 1)^ (6 + a + n + 2) (26 + 2a + 2 n - 1) 
(26 + 2a + 2n + 1)^ (26 + 2a + 2 n + 3) , for f3 = 1 

3(6 + a + 4n - 4)(6 + a + 4n - 3)(6 + a + 4n - 2)^ 
(6 + a + 4n - 1)2(6 + a + 4n)(6 + a + 4n + 1), for (3 = A 

r + for /3 = 1 
\ - for P = A ■ 



3q (g + 1) (g - 3) (^ + 4 ± 4y^TT) 



Proof. For instance, in the Jacobi case, one computes 



-'n+2 



j(4) 
-'n+1 

j(4) 



fn+3 



+ a + 6)r(^ + a + 6)r(^ + a)r (^ + 6) 
r(n + a + 6 + |)r(n + a + 6 + 2) 



r(^ + a)r(^ + 6)r(^)r(^) 

r(n + a + 6 + |)r(n + a + 6 + 3) 

^^4n+a+b r(2n + a + 6)r(2n + a + l)r(2n + 6 + l)r(2n + 3) 
r(4n + a + 6)r(4n + a + 6 + 2) 



' \/gTT = 2n + 26 + 2a + 1 for /3 = 1 and VgTT = 4n + fe + a-lfor/3 = 4 
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(„l)2 j(P) j(/3) 



(n-2)!(n + 2)! (j(/3))2 



n! 



n2 



r(/3) r(/3) 
-'n-l-'n+l 



(„_l)!(n+l)! (j^/3))2 



/3=1 



/3=4 



Q_ 



Q_ 



Qt 



13=1 



/3=4 



D 



4. Proof of Theorems 0.1, 0.2, 0.3 

From Theorems 2.1 and 2.2, the integrals In{t,c) , depending on /? = 
2,1,4, on t = (ti,t25---) and on the boundary points c = (ci,...,C2r) of E, 
relate to r- functions, as follows: 

» n 

(4.0.1) In{t,c) = / |A„(z)|^n(e^^*''^/^(^fc)^^^' 

n!rn(t,c), n arbitrary, (3 = 2 

n\Tn{t,c), n even, (3=1 

'n'lT2n{t/2,c), n arbitrary, (3 = 4. 

In{t) refers to the integral (4.QJ,) over the full range. It also follows that Tn{t, c) 

satisfies the KP-like equationll3 

(4.0.2) 



12 



Tn-2{t,c)Tn+2{t,c) p 



where 



rn{t,c) 



(KP),F : = 



5^ = (KP)t log r„(t,c) 



n arbitrary for (3 = 2 
n even for /? = 1,4 






4.1. /3 = 2, 1. Evaluating the left-hand side of (4.0.2)(for (3 = 1) yields, 
taking into account P„ := Pn{E) = In{0,c)/In{0)'- 



12 



Tn-2{t,c)Tn+2{t,c) 



Tn{t,c) 



t=0 



12 



12 



{n\f In-2{t,c)In+2{t,c) 



t=0 



(n-2)!(n + 2)! In{t,cf 



n{n - 1) I„_2(0)/„+2(0) Pn-2Pn+2 



(n + l)(n + 2) /„(0)2 

, r,, (1) Pn-2{E)Pn+2{E) 

Pl{E) ' 



p2 



l^Remember ^f ^ = 1 for /3 = 1, 4, and = for /3 = 2. 
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with hn given by Lemma 3.1. Concerning the right-hand side of (4.0.2), 
it fohows from Section 2.1 that F„(t;c) = log I„(t;c), as in (4.0.1), satisfies 
Virasoro constraints, corresponding precisely to the situation of Sections 3.1 
and 3.2 for Gauss, Laguerre and Jacobi. As explained in (3.1.4), (3.2.1) and 
(3.2.2), we express 



d^F 



dt\ 



d^F 



t=0 ^^2 



d^F 



t=0 



' dtidts 



d^F 






F = logInit,c), 



t=0 



in terms ofDk and Bk, which when substituted in the right-hand side of (4.0.2), 
i.e., in the KP-expressions, leads to (upon comparing the expressions (1.3.4) 
and (1.3.8) with (3.1.1) for Gauss and Laguerre and using (3.2.1) directly for 
Jacobi): 

{71,-1 = -f) 71,0 = 71 = 0, (5i = -§ 
72,-1 = 0,72,0 = -1/2,72,1 = 0,72 = -fcJi,(52 = 
73,-1 = -iO-1,73,0 = 0,73,1 = -^,73,2 = 73 = 0,(53 = -jCTi. 
(KP),logr„(t,c)|,^o 

= {Vj + 6nVJ + 32?! - 3P2 - ^ViV^)F + 6{VJFf + j{2 - P)n{n - 1) 
= ^((^-1 + 8(ri + (2 - /3)(n - l))Bl^ + 12^2 ^ 34^0 - l6B-iBi)F 

+6{B^^Ff + 12(2 - /3)n{n - 1] 



Laguerre with < 



7i,-i = 0,71,0 
72,-1 = 0,72,0 
73,-1 = 0,73,0 

L 73,2 = -1,73 = 



-1,71 = -f(o-i + a), 

-0-1,72,1 = -1,72 = -to-i(o-i + a), 

-Cri(T2,73,l = -(T2, 
-§CJiCr2(o-i +a). 



(KP),logr„(t,c)|^^o 



Vf - 2(/3 - 3)Vl 



2n{n - l)(/3 - 2)(/3 - 1) + (/? - 2)(4an + 4n + 5) - 4n^ - 4an - 1) P 



3 (/? - 2) /3n^ - /3n + 2an + 2n + 1 Pi 



+3Vl + 6(/3(n - 1) + a + 2)^2 - 6P3 - 4PiP3 F^ 
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-3(/3 - 2)(PiF„)2 + 6(p2 logrjv)' - 4(/3 - 3)(PiFn)(pfF, 
--(/? - 2)n(n - l)(/3n - 2/3 + 2a + 2)(/?n - /? + 2a + 2) 



- ((/3 - 2) (3 (/? - 1) (n - 1)2 + 3n2 + 6an - 4a + 2) + (a^ - l)) B^_^ 
+3 (/3 - 2) ((/3 - 1) (n - 1)2 + n^ + 2an - a) S_i - ABiB^i - 2Bi 

+2 (/3n + a) Sq-S-i + 2,BI - {pn + a) Bq] F 



+6{Bl^Ff + 4 (/? - 3) {B^iF){Bl^F) + 3 (2 - /?) {B^^Ff 
--{(5 - 2)n{n - l){fin - 2(3 + 2a + 2){f3n - (3 + 2a + 2). 
• Jacobil^ 

for 13 = 2, 
i(Z(g2-4)(KP),logT„(t,c)|,^o 

2Bt^ + {q-r + A)Bl^ - {AB-iF - s)B-i + 3gi3^ - 2qBo + SBqBI^ 

4{q - l)BiB-i + (AB-iF - s)Bi + 2(4i3_iF - s)^oS-i + 2*7^32 V 
+4Bl^F (2B0F + 3^2 ^f) 



for P = l, 
Q^ (KP), log T„(t,c) 1,^0 

= ((7 + 1) ("49^11 + 12(4e_iF - s)B^_^ + 2{q+ 12) (4^_iF - s)^oS_i 
+ 3g2^2 _ 4 (g _ 4) ^^^^_^ + g(4^„iF - s)Bi + 20(7^0^^ + 2q^B2 ]F 

+ (q2B^_, - sQiB.i + Q^B^F + ASiB.iF)'' 

- 48siB.iFf + 2Qi{B^iFf 

+ 12g2(^oF)2 + 16(7 {2q-l)B^_^FBoF + 2A{q-l)q{B^_^Ff 

+ 2i(2B-iF - s]B.iF({q + 2)BoF +{q + 3)B^_^f) + Q, 



-"^^In the Jacobi /3 = 2 case, we have bo = a — b, 61 = a + b; thus r = 2(tiQ +fci), gn = 2(2n + a + 6)^ 
and q{q^ - 4) = 16(2n + 7 + 5)2(2n + 7 + 5 - l)(2n + 7 + <5 + 1). 
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where the Qi^Q2,Q^^Qa,Q are given by (0.3.1) and where the auxUiary Qg 
happens to be exactly the one of Lemma 3.1. This estabhshes Theorems 0.1, 
0.2 and 0.3 for /3 = 2, 1, at least when 6 = 1 in the exponent of the Gaussian and 
Laguerre ensembles, upon noting that S;[logP„(i?) = S^log/„(0, c)/I„(0) = 
S^logr„(0,c). 

Finally, a simple argument captures the case 6 7^ 1. Indeed, setting aE := 
[Ji [a C2i-i,a C2i] C F, for a > 0, the elementary identities 

In{t,c) = f \An{z)ffle~'^Uzk = cf \An{z)\''fle~'ldzk 

„ n „ n 

In{t,c) = / \An{z)\''l[zte-''>'dzk = C \An{z)\''l[zte~'>'dzk, 

where C{a^ b, n, f3) is a constant independent of E, lead to the same Virasoro 
constraints as in Examples 1 and 2 (§1.3), but with the following mapping for 
the differential operators 

(4.1.1) {l3^i,Bo,Bi) ^ (^^,Bo,BiVb^ (Gauss) 



(4.1.2) -^ i^B-i,bBo,b'Bij (Laguerre). 

Therefore, the equations (0.1.2) and (0.2.2) for the probabilities (0.1.1) 
and (0.2.1) are obtained by making the substitutions (4.1.1) and (4.1.2) in 
the PDEs (0.1.2)|5=i and (0.2.2)|f,=i; this process yields the precise equations 
(0.1.2) and (0.2.2), with b ^ 1. This ends the proof of Theorems 0.1, 0.2 and 
0.3 for the cases /3 = 1, 2. 

4.2. (3 = 4:, using duality. From (4.0.1), the integral for /3 = 4 is ex- 
pressible in terms of a r-function, in which t is replaced by t/2. Hence (4.0.2) 
becomes: 

,T-2n-2 (V2,c)T2n+ 2(t/2, c) 

r2„(t/2,c) 
So, the left-hand side of (4.2.1) equals (P„ := P„(S) = /„(0, c)//„(0)) 

(n!)2 In-lit, c)In+l{t,C) 



(4-2.1) 12 ^^;^ .2 = (KP)i/2(logr2„)(t/2,c). 



^^T2n~2{t/2, c)T2n+2{t/2, c) 



T2n{t/2,C) 



t=0 



12 



12 



(n-l)!(n+l)! Init,c)^ 

n In-l{0)In+liO) Pn~lPn+l 
(n + 1) I„(0)2 P„2 



t=0 



_ ■^r,t^i4) Pn-l{E)Pn+l{E) 

where bn = bn {n, a, b) is given by Lemma 3.1 and satisfies 



6W(n,a,6) = 6«(-2n,-|,-^) 
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Recall from Theorem 1.1 (1.1.4) that In (t,c; ai,bi) and /A {t,c;ai, hi) (where 
we indicate the explicit dependence on the coefficients Oj and hi of p'/p) satisfy 
the same equations, with altered parameters: 

^k - vi^l{t-n,aiA)) li^\t,c;a^,bi) = 0, 



(^B,-Vi^l{-^t;-^n,a,,-^k)yWP\t,c;a,A) 



0. 



Setting /3 = 1 in the equations above, extracting i-partials in terms of 
BkS, and using the procedure explained in this section, we have that 

(KP),(log/«(t,c;a„6i))_^ = i?(S;n, a^, 5,) log/«(0, c; a„ 6,) 



i=0 



R{B; n, ai, h,) log pIi^\e), 



t=o 



(KP),/2(log/W(t,c;a„6,))^^g = {KP)_,^^{logli^\t,c;a„h,)) 

= RiB;-2n,ai,-b,/2)logli^\0,c;a^,h^) 
= R{B--2n,ai,-h/2)\ogP^*\E), 

where R{B; ai, hi,n) denotes the right-hand side of the equations (0.1.2), (0.2.2) 
and (0.3.4) for /3 = 1. The coefficients Oj and hi of the rational function —p'/p 
are as follows: the Oj and hi all vanish, except for 

Hermite cq = 1 ai = a2 = 6o = hi = 2h 

Laguerre ao = oi = 1 02 = 60 = — o hi = h 
Jacobi ao = 1 oi = 02 = —1 ho = a — h hi = a + h; 

thus the map 

{n,ai,bi) — > {-2n,ai,-hi/2) 

translates into the map 

(4.2.2) (n, a, b) — > {-2n, -a/2, -6/2), 

which shows that the PDEs (0.1.2), (0.2.2) and (0.3.4) for the case /3 = 4 are 
obtained by means of the map (4.2.2) from the same PDEs for /3 = 1. But 
according to (0.0.5), this is the precise way the coefficients q, s, Q^i, Qo, Qi, Q2, 
Qs; <34; Q; cvaluatcd at /3 = 4, are obtained from the same coefficients at /? = 1. 
This ends the proof of Theorem 0.3. D 

4.3. Reduction to Chazy and Painleve equations (/? = 2). Setting E = 
[-oo,x],E = [0,x],E = [-l,x] in the PDEs (0.1.2), (0.2.2) and (0.3.4) re- 
spectively, leads to the equations (0.4.1), (0.4.2) and (0.4.3) respectively, as 
announced in Section 0.4. Furthermore setting (3 = 2, the inductive terms on 
the left-hand side of (0.4.1) and (0.4.2) vanish and one obtains the ODEs: 
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• Gauss: P„(maxAj < x) = exp(— / f{u)du), where / satisfies: 

« Jx 

f" + 6 /'2 + 4fe(2n - bx^)f' + Ab'^x / = 0. 

• Laguerre: P„(maxAj < x) = exp — / du I, where / satisfies 

i \ Jx U J 

x^f" + xf" + 6x/'2 - 4//' - ((a - bxf - Anbx)/' - b{2n + a- bx)f = 0. 

• Jacobi: P„(max Aj < x) = exp ( — / Trdu ] , where / satisfies: 

i \ Jx l-U^ J 

2(x2 - iff + 4(x2 - 1) {xf - 3/2) + (l6x/ - qn{x^ - 1) - 2sx - r) f 

-f (4/ - qnX -s)=0, 

where r, s, Qn are defined in (0.3.1). 

These three equations are of the form 

(4.3.1) /"' + 5/" + l^r - '-^ir + J/^ + ^/' - ^/ + ^ = 0. 

with the following coefficients P,Q,R: 

Gauss P(x) = 1 4Q(x) = -Ab^x^ + 8bn R = 

Laguerre P{x) = x 4(5(x) = — (fex — af + 46nx R = 

Jacobi P(x) = l-x2 4Q(x) = -^(^^(x^ - 1) + 2sx + r) i? = 0. 
The general Chazy class of differential equations are equations of the form 

/ = F{z, f, f , f ), where F is rational in /, /', /" and locally analytic in z, 

subjected to the requirement that the general solution be free of movable 
branch points; the latter is a branch point whose location depends on the 
integration constants. In his classification, Chazy found thirteen cases, the 
first of which is given by (4.3.1), with arbitrary polynomials P{z),Q{z),R{z) 
of degree 3,2,1 respectively. 

Cosgrove ([11], [12]), (A. 3), shows this third-order equation has a first 
integral, which is second-order in / and quadratic in /", 

(4.3.2) f +^ [{Pf + Qf + R)r - {P'f + Q'f + R')f 

+liP"f' + Q")f-lp"'f + c)=o, 
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with an integration constant c. In the three cases, discussed above, c = 0. 
Notice equations of the general form 

f"^ = Gix,f,f') 

are invariant under the map 

aiz + 02 , „ as/ + qqz + ay 

X 1-^ and / 1-^ . 

a^z + a4 a^z + a4 

Using this map, the polynomial P{z) can be normalized to 

P{z) = z{z — 1), z, or 1. 

Equation (4.3.2) is a master Painleve equation, containing the six Painleve 
equations. If /(x) satisfies the first three equations above, then the new func- 
tion g{z), defined below. 

Gauss g{z) = h'^l^ f{zh'^/^) + fnz 
Laguerre g{z) = f{z) + \{2n + a)z + ^ 



Jacobi g{z) := -kf{x)\x=2z-i - iz + 



q+s 
16 



satisfies the following canonical equations of Cosgrove and Scoufis ([11], [12]): 
• /2 ^ _4^/3 _^ ^^g> _ g^2 j^ ^^g> j^ ^^^ (PaiulevB IV) 

. {zgy = {zg' - g) {-Ag'^ + A^{zg' - g) + A^ + A.J + A^, 

(Painleve V) 

. {z{z - l)g"f = {zg' - g) [^g'^ - Ag'{zg' - g) + A^ + A^g'^ + A^g' + A^ 

(Painleve VI) 
with respective coefficients 

• ^l=3(f)^^2=-(f)^ 

. Ai = h\ A2 = }?{{n + 1)2 + 4), ^3 = -a%{n + f ), ^4 = ^ 

•((n + f)' + ^), 

. Ai = ^, ^2 = f , ^3 = ^2^^^, A, = ^(2.2 + ,.). 

Each of the equations above can be transformed into the standard Painleve 
equations. 
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5. Appendix. Self-similarity proof of 
the Virasoro constraints (Theorem 1.1) 

Given the data (0.0.1) to (0.0.3), namely p = e~^ and —p'/p = V = 
alf = EcT biz'/Y.'^ aiz' and ^ = Ui [c2i-i,C2i\ C F C M, we show that the 
muhiple integral 

„ n 

(5.0.1) /„(t,c;/3) := / \l^n{x)f \{ (eEr*«-lp(xfc)dxfe) , forn > 

satisfies the Virasoro constraints of Theorem 1.1, using a (much less concep- 
tual!) self-similarity argument. Setting 

n 

dln{x) := \Anix)f n (e^^*'"'^p(xfc)d: 



fc=i 



IXk , 



we state the following lemma: 



Lemma 5.1. The following variational formula holds: 



(5.0.2; 
d_ 
de 



dln{xi ^ Xi + e/(xi)xf+-^) 



e=0 



E («^ ''4'i,n - h ''j£,+l,„) din. 



i=0 



Proof. Upon setting 

n 

(5.0.3) E{x,t) := l[e^ti'^^lp{xk) 



-Q g-V(xfe,t) ^ ^j^gj^g y^^^ ^^ ._ y^^-^ -Y,t, 



X 



the following two relations hold: 
(5.0.4) 

\ 

'Jk,o 




--5 
dhdti 2 ''' 



J 



J2 ^-^^« + 



k-i 



'^a-^P 



l<a<l3<n 
\ i,3>0 



l<a<n 



E, 



W, + "^'^'^ ' ^ 



J2 x^^] E, for ah k>0. 



. l<a<n 

So, the point now is to compute the e-derivative 

(5.0.5) -^(|A„(x)|^eELi(-V'(-.)+E.^i*»-l)dxi...dx„, 

which consists of three contributions: 



e=0 
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Contribution 1: 



(5.0.6) 
d 



de 



^{x + ef{x)x^+') 



e=0 



/3|A(x)|^ E -7r^og(\x^-x^ + e{f{x^)xi+^-f{x^)x: 



l<a<7<n 



de 



fc + lN 

7 ' 



e=0 






l<a<7<n 



/3|A(x)|/5^o, ^ 



ol ■7 



^=0 l<Q<7<n ^" ^T 



/ 



/?|A(x)|'^E«^ 
£=0 



E a;'„x^ + (n - 1) ^ ^^^^ 5 '^^+'^.0 



\ i,j>0 

\ l<a<7<n 



l<a<n 



/3^-i|A(x)|/^Ea4i Y. ^-f5.+.,o 

i,j>0 

( k+l+l\ ( d \ n(n-l) \ 

00 
/3i?-i|A(x)|^^a, 

A V- ^^ / A: + £+l\ a n(n-l), \^ 



2 ,^^^, 9t,9t, 

i,3>0 



dtk^ 



Contribution 2: 
(5.0.7) 



nc^(a;a + e/(Xa)x^+^) 



ae 



£ = 



E (/'(^«)^^' + (A;+ i)/(x«)x^) n^^ 



1 

00 



n n 



1=0 a=l 1 

d 



£=0 



ae. 



at 



k+e 



+ n6k+e,o I eJJ dxi 
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on / 

^ a=l V 



Contribution 3: 

oo 
a=l i=l a=l / 



oo n 

i=l a=l 

n oo n 



e=0 



oo n n 

E^^E^^'^-' + E^^^^^E 

i>l 



^i+k+e 



d 



T.be(j^-— + n5, 



— ^9tk+e+i 

oo oo 



k+e+1,0 



+ E '^^ E ^*« ( ^I ^ n6i+k+i,o ] E. 

i=0 ^=l \Oti+kU J) 

As mentioned, to conclude (5.0.2), we must add up the three contributions 
(5.0.6), (5.0.7) and (5.0.8), resulting in: 

(5.0.9) 
d 

— dln{xi ^ Xi+ £f{Xi)x^^'^] 

OS £=0 



(e«^ (^j5 + M + (^ + A: + i)(i - f ))4'i 



+n{{n - 1)- + l)5k+ifi\ - E ^^ (4+^+1 + nSk+i+ifij dln{x). 

where j], := ^J^, as in (1.1.8). Thus we use (1.1.8) to end the proof of 
Lemma 5.1. D 

Proof of Theorem 1.1. The change of integration variable Xi i-^ Xi + 
ef{xi)x^~^ in the integral (5.0.1) leaves the integral invariant, but it induces a 
change of limits of integration, given by the inverse of the map above; namely 
the Cj's in i? = Ui[c2i-i,C2i], get mapped as follows: 

Q^c,-e/(Q)cf+^ + 0(e2). 
Therefore, setting 

E^ = U[c2.„i - e/(c2._i)4t\ + 0(e2),c2. - e/(c2.)4+' + 0(8^)], 
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we find, using Lemma 5.1 and the fundamental theorem of calculus, 



2r 



This ends the alternative proof of Theorem 1.1. D 
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